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Motivated by recent experiments on quantum trapped fields, we give a rigorous proof tliat repeated 
indirect quantum non-demolition (QND) measurements converge to tlie collapse of the wave function 
as predicted by the postulates of quantum mechanics for direct measurements. We also relate the 
rate of convergence towards the collapsed wave function to the relative entropy of each indirect 
measurement, a result which makes contact with information theory. 
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Wave function collapse is a basic axiom of quantum 
direct measurement a la Von Neumann [i]. A quantum 
non-demolition measurement [5] is one for which the col- 
lapsed state is an eigenstate of the free evolution. Repeat- 
ing the measurement on the collapsed state yields identi- 
cal results since this state is preserved by the evolution. 
Indirect measurements [6j consists in letting the quantum 
system under study be entangled with another quantum 
system, called the probe, and in implementing a direct 
measurement on the probe. Since the system and the 
probe are entangled, one gains information. Repeating 
the process of entanglement and measurement increases 
statistically the information one gets on the system. 

Developing experimental and theoretical expertise on 
quantum measurement processes is mandatory for devel- 
oping quantum state manipulation. It was early realized 
[3 [5] that modeling quantum measurements require sys- 
tems with infinitely many degrees of freedom, e.g. as in 
the phenomenological stochastic models of ^ . The need 
to describe quantum jumps and randomness inherent to 
repeated measurements lead to the concept of quantum 
trajectories jTUHH]. In parallel, tools of open quantum 
systems, specifically those of quantum stochastic calcu- 
lus [H] , have been adapted to the description of quantum 
continual measurements [13] and quantum feedback [13] . 
In most of these stochastic models, the driving noises, 
often classical or quantum Brownian motions, are linked 
to the degrees of freedom of the measurement apparatus. 
Although bearing similarities with these frameworks, our 
proof of the wave collapse in series of QND measurements 
is based on a purely quantum description of the repeated 
probe-system interactions. 

Experiments on repeated indirect quantum non- 
demolition measurements have recently been performed, 
in particular in quantum optics. As an example, let us 
look at [TS] whose setup is the following. The tested 
quantum system is a resonant electromagnetic cavity se- 
lecting photons of given frequency. It is probed by send- 
ing Rydberg atoms through it, one after the other. Dur- 
ing the atom-photon interaction each atom behaves as a 



two-state system modeled by a spin one-half [16 . The 
atoms are prepared with their effective spins pointing in 
the Ox direction |17| . The experimental protocol ensures 
that the atom effective spin rotates around the Oz axis by 
an angle proportional to the number of photons nph in 
the cavity, say riph^ with 9 a fixed angle. After interac- 
tion, the atom-photon system is entangled, but the cavity 
state gets unchanged if it is initially an eigenstate of the 
free photon hamiltonian. The effective atom spin is then 
measured along a direction perpendicular to Oz but at an- 
gle (f) with respect to Ox. The output of the spin measure 
is ± with probabilities p+(0|nph) = cos^[(nph^ — ^!>)/2] 
and p_(0|nph) = sin^[(nph6' — 0)/2], if there are nph 
photons in the cavity. If the initial photon distribu- 
tion is qoifiph), the probability to measure an effective 
spin ± is J2h h 9o('^ph)p± (0|"-ph)- No direct measure- 
ment on the cavity is done. The experimental aim is 
to reconstruct the initial photon distribution by accu- 
mulating informations from the repeated atom effective 
spin measurements. The photon distribution is recalcu- 
lated after each atom measurement using Bayes law |18j . 
Figjl] shows experimental data for the evolution of re- 
constructed photon distributions. For each realization, 
they converge, experimentally and numerically |15j . to 
peaked distributions whose centers depend on the real- 
ization. This is the collapse. 

Let us abstract and generalize the previous situation. 
At initial time, the system is in state |(^o) = Iv)- It inter- 
acts during time At with a probe initially in state \4>), so 
that the pair (probe-f system) evolves into U{\4>) <E) |<p)), 
where U is some unitary operator acting on the Hilbert 
space Hprobe "8) Hsyst- After At, the system-probe inter- 
action can be neglected. A perfect measurement a la Von 
Neumann is then performed on the probe. This means 
that there is an orthonormal basis i G /, of Uprobe 
such that, after the measurement, the (probe-|-system)- 
state is proportional to (|z)(z|(8)Id)C/(|'i/')<8'|v)) with prob- 
ability ||(|i)(i|(g)Id)t/(|^/')(g)|(p))||2. The vanishing of this 
probability for a certain state \i) means that the probe 
cannot be found in state \i), so we can (and shall) simply 
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FIG. 1: Two experimental samples ]19l of reconstructed pho- 
ton distributions as functions of the number of indirected mea- 
surements (i.e. the number of atoms traversing the cavity) 
according to \15^ . The collapse of the photon distribution to a 
realization dependent sharply defined number is clearly visible. 

forget about that possibility. 

We make the fohowing assumption, related to non- 
demolition, on the evolution operator U: there is an or- 
thonormal basis \a), a A, of Hsyst and a collection of 
operators Ua acting on Uprobe such that, for each a, 

U{\t^)(g>\a)) = {U^\^))(^\a). (1) 

The operators Ua are automatically unitary. If the probe 
is found in state \i) after the measurement, (probabil- 
ity EqsA l(*|f^"l^)Pl("IV'>P)> the pair (probe-Hsystem) 
is again in a tensor product state \i) (g) \(pi) where 

l*^!/ = ~, ZT72- ^'^^ 

{EaeAmUaimi^MT 

It is clear that the motivating experiment fulfills this 
property if |a) is the occupation number basis |20| . 

The physics of this hypothesis is that the final aim is 
to measure an observable on the system for whom the 
states \a) are eigenstates. As we shall see, this (direct) 
measurement can be (indirectly) achieved by repeated 
measurements on successive probes. So, one presents 
another probe to the system in state \(pi), let them 
interact and, after interaction, measures the probe to 
get \(f2) and so on. Notice that, in general, at each step 
one could change the probe initial state, the observable 
measured on the probe (this is indeed what happens in 
the motivating example), and even the type of probes: 
the only thing one has to keep fixed is the basis \a) 
for which property ([T]) holds. Most of the following 
discussion can be extended to the general setting [21] 
but to keep notation simple, we concentrate on the case 
when 1-0) and the basis \i) are the same for all probes. 

We start with a summary of our results: 
i) If a series of repeated indirect measurements is 
conducted, the state of the system will stabilize over 
time and go to a limit. Carrying identical independent 
experiments again, the system state will stabilize over 



time again but possibly with different limits. 
a) Under a physically meaningful non-degeneracy 
condition, the only possible limits for the state of the 
system are the pointer states |a), and the probability to 
end in state |a) starting from state is Hence 
the outcome of a large number of repeated indirect mea- 
surements satisfying condition ([T]) obeys the standard 
rules of quantum mechanics direct measurements. 
Hi) Under the same non-degeneracy condition, the 
measurements on the probes allow to infer the limit 
pointer state for each independent experiment. 
iv) The rate of convergence to one of the pointer 
states is governed by the relative entropy of certain 
probability measures in classical probe space. The order 
of magnitude of the probability that, while the repeated 
measurements are conducted, the state of the system 
comes close to a pointer state but ends up finally in 
another one can be computed explicitly. 

The tools to prove these statements come from the 
classical theory of random processes : strong law of large 
numbers, martingale convergence theorem, large devia- 
tions. A proof of the wave function collapse using the 
martingale convergence theorem appeared in [22] . These 
works are based on non-linear stochastic extensions of the 
Schrodinger equation |23j whereas our results are pure 
consequences of quantum mechanics (with measurements 
on probes) [24] and are closer in spirit to quantum tra- 
jectory approaches |10l 111) and to experiments. 

We now turn to the proofs. One can rephrase eq.([2]) 
by saying that, for each a € A, 

= im.Wi^l^o) 

if the probe is found in state Thus, a crucial conse- 
quence of ([I]) is that there are no interference terms for 
different a's, so that taking the modulus squared does 
not lead to (much) loss of information. We set p{i\a) = 
mU^m', and qo{a) = \{a\ipo)\^ qi{a) = \{a\ipi)\^ 
(72(a) = |(a|(/J2)P and so on. Observe that after measur- 
ing the n'^ probe one has, for each a £ A, 

g„+i a) = g„ a = (3) 

with probability 7r„(i) = J2i3(^A9n{l3)p{i\l3). 

This is a random recursion relation which is of Marko- 
vian type : to compute the possible values of (?„+i(a) and 
their respective probabilities, all one needs to know are 
the (?„(/3)'s. Each probe measurement leads to a choice 
among the probe states \i) such that 7r„(i) ^ 0. The 
question to be settled is the long time behavior of the 
resulting random sequences (?„(a). 

Observe that the (7„(q!)'s and the p(i|a)'s are > 0. 
Moreover, J^a^A 9n(a) = 1 and I]ie/^'(*l") = ^ ^^^^ 
a G A. It follows that X^ie/ = 1 as it should be. 
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A crucial question is the following : having observed the 
random sequences gm(/3) for m — 0, - ■ ■ ,n and all /3's in 
A, what is the average value of qn+i{a)7 From ([3|, it 
is immediate that this (conditional) average, which we 
denote by E(g„+i(a)|(7o, • • • , ^n), is 

E{qn+i[a)\qo,--- ,qn) = 2^ q„[a) 7TT-7r„(i) 

= X! qn{a)p{i\a). 

i,7r„(i)#0 

Now, 7r„(i) = X^^eA ^nd^)?'!*!/^) ^'^'^ '^'-"^ ^^^^ vanish, 
the product qn{0)pii\P) has to vanish for all P & A, and 
in particular for /3 — a, so that ^ {i)^o 1nict)pi'i\c() = 
X^ig/ — Qnic()- Hence we find that 



E((7„+i(a)|go, ■ • • = quia). 



(4) 



In the theory of random processes, such a property 
defines the concept of martingale : the sequence qo,qi, - ■ ■ 
is a martingale, because if one knows it up to time n (i.e. 
if one knows qo, ■ ■ ■ ,qn) its expectation at time n+1 is its 
value at time n (i.e. g„). To connect quantum measures 
to conditional expectations is not so surprising because 
both rely on orthogonal projections in Hilbert spaces. 

The martingale at hand has a peculiar property: it is 
bounded (every (?„(a) is > and J2a€Alni'^) ~ We 
can then quote a special case of the martingale conver- 
gence theorem (see any modern textbook on probability 
theory, e.g |27j , for a precise statement): A random se- 
quence <7o: 917 ■ ■ ■ which is a bounded martingale converges 
almost surely and in L^. The limit, a random variable 
Qoo, is such that its expectation satisfies E{qrx,) — qo- 

This is a deep theorem and there is no intuitive argu- 
ment that we know to explain it |25j . But in our case 
its meaning is simple. The statement of almost sure con- 
vergence is precisely the mathematical formulation of i). 
The statement of convergence is a simple consequence 
of the Lebesgue dominated convergence theorem, because 
our martingale is bounded. The statement on the expec- 
tation of the limit random variable yields the second part 
of ii) once we have given an independent argument to 
show that the possible limits are the pointer states. 

To get this, we observe that the convergence of qn{a) 
leads to the convergence of 7r„(i) = J2p<£A^n{P)p{i\f3). 
If i is such that 7roo(j) 7^ then, for n large enough, 
7rn(*) > 7roo(i)/2 > which implies that, with probability 
1, the n*'' probe will be found in state i for arbitrarily 
large values of n. This allows to take the large n limit in 
([3| for this value of i. Hence 



goo (a) = goo (a) 



p{i\a) 



E/3eA9oo(/3)p(i|/3)^ 



to these a's. Then, we can simplify to get p{i\a) = 7roo(i) 
for any i such that 7roo(«) ^ 0. If goo (a) 7^ 0, 7roo(j) = 
implies p{i\a) = 0, so that p{i\a) — 7roo(i) is actually 
valid for any i. The right-hand side may depend on i but 
it does not depend on a. So the same holds for the left- 
hand side: this means that the evolution operator U and 
the probe measurement act in a degenerate way on the 
corresponding kets |a). In such a degenerate situation, 
we cannot expect to measure them individually, just as in 
a standard quantum measure of a system observable we 
cannot separate the |q;)'s having the same eigenvalue |26j . 
So, we assume that for any a, (3 E A there is some i G I 
such that p{i\a) ^ p(i\/3), and we get that goo(a) = ^a.-y 
for some 7, i.e. the only possible values for goo(Q;) are 
or 1. The equality E{q^{a)) — go (a) then implies that 
goo(Q;) takes value 1 with probability go (a) — \{a\ip)\'^ 
and with probability 1 — go (a) as expected in a perfect 
measurement of a non-degenerate system observable with 
the |a)'s as eigenstates. 

The proofs of statements Hi) and iv) use the same 
tools. We start by determining the rate of convergence 
to the limiting system state. This turns out to depend on 
this limiting state and this is also the clue to statement 
Hi). What we have proved so far implies that at some 
time, say no, one of the components, say gno(7): "^iH be 
large, i.e. close to 1, so that all other components will be 
small. We can then replace ([3| by an approximate linear 
recursion relation, namely, for a ^ j, 



qn+iice) = qn{a) 



p(i|a) 
p{i\j) 



(5) 



with probability p{i\"f) (if non zero). The proof given 
above shows again that this random recursion relation 
defines a martingale. There is a subtle point however : 
this martingale is not bounded anymore and the martin- 
gale convergence theorem does not apply. However, we 
can rely on a simpler tool. Defining In = logg„ we get, 
for a 7^ 7, 



In+iia) = ln{a) + lo| 



,p{i\a) 
P{i\l) ' 



(6) 



for any i G I such that 7roo(j) 7^ 0. Only the a's for which 
goo (a) 7^ yield a nontrivial equation, so we can restrict 



with probability p(i|7) (if non zero). So Z„(a) — 
Zo(a) is the sum of n independent identically dis- 
tributed random variables with mean —S^jla) = 
XjP(i|7) logp(i|Q;)/p(i|7). Remember that for each /3, 
the collection p{i\f3), i G I, defines a probability on /, 
and S'(7|q:) is nothing but the relative entropy of p{i\j) 
with respect to p{i\a), a quantity which is always non- 
negative, and in fact strictly positive under the non- 
degeneracy assumption. The law of large numbers yields 
ln{ct) ~ —nS{'y\a) —00, so that g„(a) converges expo- 
nentially to with rate S'(7|q;). Hence, as soon as one of 
the components, say g(7), has become reasonably close 
to one, with high probability the state of the system will 
converge to I7). In this situation, each measurement on 
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the probe leads to a gain of information on the system 
state which in average is given for each component a 7^ 7 
by the relative entropy S{^\a). 

By the strong law of large numbers, the previous dis- 
cussion also implies that if the limit state is I7), the fre- 
quency of measurements leading to probe state |i) will 
converge to p{i\j). By the non-degeneracy hypothesis 
this fixes the limit pointer state unambiguously. This 
proves statement in). In practice, an histogram of all 
rii/n, the fraction of probes measured in state \i) in a sin- 
gle series of a large number n of repeated measurements, 
for i e /, will be close to p(i|7) for a single I7) , allowing to 
identify I7). Then conducting many independent homo- 
geneous series (starting each experiment with the same 
system state) allows to reconstruct the probabilities go. 
Hence the homogeneous repeated measurement scheme is 
fully equivalent to an ideal Von Neumann measurement. 

To finish the discussion, note that by the martingale 
property, knowing the results of probe measurements up 
to time no I the probability to end in pointer state I7) 
is exactly 9no(7), which is close to 1. The quantity 
1 — (7) is the probability to end in another pointer 
state. It is also the order of magnitude of the probability 
that the above discussion breaks down. This occurs pre- 
cisely when the random evolution invalidates the linear 
approximation. If this happens, it will be likely to hap- 
pen quickly after uq because, if for a long time after hq 
the g„(a) remain small, the law of large numbers implies 
that they are very likely to decrease exponentially so that 
escaping away from the pointer state I7) will get harder 
and harder. Take some e > such that if 1 — e < 9,1(7) 
the linear approximation is good to describe the transi- 
tion from time n to time n -|- 1. Suppose that during a 
random evolution this condition on qn{l) remains valid 
for riQ < n < ni. By standard large deviation theory 
(Cramer's theorem) if rii — uq is large, the probability 
that, for a given a, {a) is of order e (instead of being 
of order £exp[— (ni — ?io)S'(7|a)]) is estimated crudely 
as ^ A"^""'' for a certain A* < 1 which is the minimum 

over s > of the function A(s) = Y^iViAl) (f^) ■ 

Finally, we emphasize that the (infinite) series of 
indirect experiments may be viewed as building a 
measurement apparatus [21] • Indeed, the reading of 
the asymptotic behavior of the frequencies of the probe 
measurement outcomes allows to register the limit 
pointer state. 
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